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Abstract 

We describe the formalism, and present the results, for a triple-Regge analysis of the 
available pp and pp high-energy data which explicitly accounts for absorptive corrections. 
In particular, we allow for the gap survival probability, S 2 , in single proton diffractive 
dissociation. Since for pp scattering the value of S 2 is rather small, the triple-Pomeron 
^ ■ vertex obtained in this analysis is larger than that obtained in the old analyses where the 

suppression caused by the absorptive corrections was implicitly included in an effective 
vertex. We show that the bare triple-Pomeron coupling that we extract from the pp and 
pp data is consistent with that obtained in a description of the jp — * J /if) + Y HERA 
data. The analyses of the data prefer a zero slope, corresponding to the small size of the 
bare vertex, giving the hope of a smooth matching to the perturbative QCD treatment of 
the triple-Pomeron coupling. 



1 Motivation 

The total and elastic proton-proton cross sections are usually described in terms of an eikonal 
model [1]. The advantage of using an eikonal framework is that it automatically satisfies 
s-channel unitarity, which follows once we include the elastic rescattering of the interacting 



particles. Furthermore, the Good- Walker formalism [2] allows the possibility of excited proton 
states occurring in the intermediate state. In this way we include low-mass diffractive dissoci- 
ation. It was demonstrated in [3] that a two-channel eikonal is sufficient to capture the main 
features of this dissociation (provided, of course, we are not interested in the production of a 
particular resonant state), see also [4]. 

On the other hand, high-mass (M) diffractive dissociation is described in terms of the triple- 
Regge formalism, where the differential cross section, da/dtdM 2 , is driven by the triple- Regge 
couplings gR 1 R 2 R 3 . The values of these couplings were determined from the data available in 
triple- Regge domain in the 1970s; see, for example [5, 6, 7]. However, these early analyses 
did not allow for absorptive corrections, which are sizeable in hadron-hadron collisions. There- 
fore, the old triple-Regge couplings must be regarded, not as the bare vertices, but as effective 
couplings embodying absorptive effects. It was pointed out in Ref. [8] that the original bare 
triple-Pomeron vertex may be about three times larger than the effective coupling which is ob- 
tained with the neglect of absorptive corrections. However, absorptive effects are very sensitive 
to the t-slope of the triple-Reggeon vertices and in the 1970s it was "impossible to determine the 
exact values of the bare triple-Pomeron coupling" [8]; the accuracy and energy range covered by 
the data were insufficient. So the results of the triple-Regge analyses were presented in the form 
of effective couplings. This leads to a problem, since the absorptive effects are not constant 
factors, but depend on the energy and the type of diffractive process. Since the inelastic cross 
section expected at the LHC is more than twice as large as that observed at fixed target and 
CERN-ISR energies, we cannot use the results of the old triple-Regge analyses to predict the 
diffractive effects at the LHC 1 . It is therefore necessary to perform a new triple-Regge analysis 
that includes the absorptive effects explicitly. Here we present the first analysis of the avail- 
able diffractive (and elastic) data in the CERN-ISR - Tevatron region in which the absorptive 
corrections are systematically included. The absorptive corrections are calculated in terms of a 
two-channel eikonal model fitted to describe the total and elastic differential cross section data 
for pp and pp scattering. 

At the moment, the energy behaviour of the scattering amplitude may be consistently 
described by two different scenarios for the asymptotic regime [11]. One is called the weak 
coupling of the Pomerons. In this case, at very high energy, the cross sections tend to the 
universal constant value 

o"tot constant as s — > oo. (1) 

In order not to violate unitarity, the triple-Pomeron coupling must vanish with vanishing trans- 
verse momentum transferred through the Pomeron [12] 

93P oc q 2 t as q t -> 0. (2) 

Another possibility is called the strong coupling scenario [13]. Here, at a very high energies, the 
cross sections grows as 

CT to t oc (In s) v with < rj < 2, (3) 

1 Already, in order to describe data up to the Tevatron energy, Goulianos and Montanha [9] found it necessary 
to restore unitarity by renormalising the "Pomeron flux" [10]. This mimics to some extent the absorptive effects. 
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and the bare vertex 



constant. 



(4) 



The present data are usually described within the Froissart-like limit of the second scenario 
(with rj — 2). However to reach asymptotics we need very high energy - the energy at which 
the slope of the elastic amplitude, B = B Q + a' P \n(s) is dominated by the second term, that is 
when apln(s) 3> B . This is far beyond the energies available at present. Another possibility, to 
distinguish between the weak and strong approaches, is to study the q t dependence of the bare 
triple-Pomeron vertex. Thus, it is important to extract the bare vertex before its behaviour is 
affected by absorptive corrections. 

This paper is devoted to a detailed analysis of "soft" scattering data using, in turn, the 
strong and weak triple-Pomeron coupling behaviours. In Section 2 we recall the Good- Walker 
formalism for the two-channel eikonal model, in Section 3 we briefly describe the parametri- 
sation used to calculate the absorptive corrections for the triple-Reggeon cross sections and 
demonstrate that this parametrisation is consistent with the available data on the differential 
elastic (da e \/dt) and total cross sections. Next, in Section 4, we give the formulae for the 
inclusive cross section, M 2 da/dtdM 2 , in the triple-Regge region which incorporates the screen- 
ing corrections caused by the two-channel eikonal. The results of our triple-Regge analyses 
of pp — > pX and pp — > pX data are presented in Section 5. We perform fits using both the 
strong and weak triple-Pomeron coupling. We find that the data favour the former ansatz. In 
Section 6 we discuss the description of the data for inelastic diffractive J/vp photoproduction, 
7P — > J/ip + Y, obtained at HERA. We find that the J/vp HERA data again favour the strong 
triple-Pomeron coupling scenario. In Section 7 we present our conclusions. 

2 Resume of the eikonal formalism 
2.1 Single- channel eikonal model 

First, we briefly recall the relevant features of the single-channel eikonal model. That is we focus 
on elastic unitarity. Then "disc T" is simply the discontinuity of the pp scattering amplitude T 
across the two-particle s-channel cut. At high energies we have a sizeable inelastic component. 
The s-channel unitarity relation is diagonal in the impact parameter, b, basis, and may be 
written 



2ImT cl (s,&) = |T cl ( S ,6)| 2 + G incl ( S ,6) 



(5) 



with 




(6) 



(8) 



(7) 
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These equations are satisfied by 

ImT el (s,&) = 1-e"^ 2 (9) 

a el (s,b) = (1-e"^ 2 ) 2 , (10) 

a hlcl (s,b) = l-e~ n , (11) 

where Q(s,b) > is called the opacity (optical density) or eikonal 2 . From (11), we see that 
exp(— Q(s,b)) is the probability that no inelastic scattering occurs. 

The opacity corresponding to an individual Regge pole i = R or P takes the form 

Sk(s,b) = - I d\ M*,<h)e** = A(0)2(s/ ; or ' (0) " e"^, (12) 

71 J Di 

where q\ = —t, and where the contribution to the amplitude from Reggeon i exchange is 3 

A(s,q t ) = vM^is/sor^- 1 , (13) 
where the t behaviour of the Reggeon-proton vertex is taken to have the exponential form 

Pi(t) = A(0)exp(r,t/2). (14) 

Thus, for a linear trajectory 

cti = Q!j(0) + a'jt, (15) 
the t-slope of the amplitude Ai is given by 

Bi = n + a[ \n(s/s ). (16) 

For Regge exchange with signature ±1, the factor rji in eq. (13) is 

— (I ± e -™<*i) 
sm7raij 

Although we take exponential forms for the exchange of the secondary Reggeons, we choose, 
as in Ref. [14], a power-like form for the proton-Pomeron vertex 

(1 - t/ai){\ - t/a 2 ) 

Also, to improve the description of the low \t\ region, we follow Anselm and Gribov [16] and 
include pion-loop insertions 4 in Pomeron exchange which results in a non-linear form of the 
Pomeron trajectory 

a P (t) = a P (0) + *' P t-^h (^pi , (19) 



2 Sometimes Cl/2 is called the eikonal; for simplicity we omit the real part of T e \. At high energies, the ratio 
RcT e i/ImT e i is small, and can be evaluated via a dispersion relation. 

3 Note that, unlike [3, 14, 15], here we use the normalisation cr* ot = inlmrji /? i (0) 2 (s/s ) ai< ' ^ 1 , rather than 
a tot _ i mT j. ^((^(s/so) '^ -1 - ^ na ^ j S) as compared to our previous papers, here all the coupling constants 
(fii and guj) are decreased by a factor ^/4jt. 

4 Note that, in (20) below, we have corrected the misprint which occurs in the published version of [16]. 



where 



h(r) 



-Fl{t) 2r - (l + rp In 



( 



\/TT7+ 1 



v / TT7- 1 



) 
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(20) 



with r = 4to^/|£| and to = 1 GeV. The coefficient /5 2 specifies the irn total cross section, 
and F n (t) is the form factor of the pion-Pomeron vertex. The coefficient (3 2 m 2 T /32ii 3 in (19) is 
small 5 , but due to the tiny scale m n the t dependence of h{r) is steep and non-linear. It has an 
important effect on the local slope, d(\nd<j c \/dt)/dt, of the elastic cross section. For the results 
that we obtain below for the Pomeron trajectory, ap(t), it is important to note that expression 
(20) for h(r) has been renormalised [16], such that 



where h n (r) denotes the full pion-loop contribution. Further discussion of the pion-loop con- 
tribution can be found in Ref . [14] . 

Finally, the total opacity is given by the sum of all the allowed Regge exchanges 



In the present simplified analysis we neglect the real part of the amplitude, since our goal is 
to investigate the role of the absorptive effect which is caused essentially by the imaginary part 
of the amplitude, that is by the real part of the opacity. Therefore, the values of the couplings 
Pi(0), that we quote in Table 1 below, are actually Im^j/^O)), where the rji are the signature 
factors of (17). 

2.2 Inclusion of low-mass diffractive dissociation 

So much for elastic diffraction. Now we turn to inelastic diffraction, which is a consequence 
of the internal structure of hadrons. This is simplest to describe at high energies, where the 
lifetime of the fluctuations of a fast hadron is large, r ~ E/m 2 , and during these time intervals 
the corresponding Fock states can be considered as 'frozen'. Each hadronic constituent can 
undergo scattering and thus destroy the coherence of the fluctuations. As a consequence, the 
outgoing superposition of states will be different from the incident particle, and will most likely 
contain multiparticle states, so we will have inelastic, as well as elastic, diffraction. 

To discuss inelastic diffraction, it is convenient to follow Good and Walker [2], and to 
introduce states (fik which diagonalize the T matrix. Such eigenstates only undergo elastic 
scattering. Since there are no off-diagonal transitions, 



h( T ) = K(t) - K(0) 



(21) 



(22) 



i=P,R 



{<j> j \T\<f> k )=0 iorj^k, 



(23) 



5 We use the additive quark model relation /3 T (= g^p) = §/3p(0). 
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a state k cannot diffractively dissociate into a state j. We have noted that this is not, in general, 
true for hadronic states, which are not eigenstates of the S'-matrix, that is of T. To account for 
the internal structure of the hadronic states, we have to enlarge the set of intermediate states, 
from just the single elastic channel, and to introduce a multichannel eikonal. We will consider 
such an example below, but first let us express the cross section in terms of the probability 
amplitudes F k of the hadronic process proceeding via the various diffractive eigenstates 6 <p k . 

Let us denote the orthogonal matrix which diagonalizes ImT by a, so that 

ImT = aFa T with {4>j\F\4> k ) = F k 5 jk . (24) 
Now consider the diffractive dissociation of an arbitrary incoming state 

\j) = E l&>- (25) 

k 

The elastic scattering amplitude for this state satisfies 

(j\lmT\j) = £ \a jk \ 2 F k = (F), (26) 

k 

where F k = (<f) k \F\(f) k ) and where the brackets of (F) mean that we take the average of F over 
the initial probability distribution of diffractive eigenstates. After the diffractive scattering 
described by Tfj, the final state |/) will, in general, be a different superposition of eigenstates 
from that of \j), which was shown in (25). At high energies we may neglect the real parts of 
the diffractive amplitudes. Then, for cross sections at a given impact parameter b, we have 



da. 



21m(j\T\j) = \a jk \ 2 F k = 2(F) 



d?b 

d(J e l 



\(J\T\J}\ 2 = £ \a jk \ 2 F k ) = (F) 2 (27) 



d 2 b 

= E mm 2 = E M 2 ^ 2 = {f 2 ). 
a k k 

It follows that the cross section for the single diffractive dissociation of a proton, 

(F 2 ) - (F) 2 , (28) 



d<7 SB 



d 2 b 

is given by the statistical dispersion in the absorption probabilities of the diffractive eigenstates. 
Here the average is taken over the components k of the incoming proton which dissociates. If 
the averages are taken over the components of both of the incoming particles, then in (28) we 
must introduce a second index on F, that is F ik , and sum over k and i. In this case the sum is 
the cross section for single and double dissociation. 



6 The exponent exp(— £l k ) describes the probability that the diffractive eigenstate (f> k is not absorbed in the 
interaction. Later we will see that the rapidity gap survival factors, S 2 , can be described in terms of such 
eikonal exponents. 
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3 Description of the elastic data 



In preparation for the triple- Regge analysis of the data for high-mass diffraction, we first perform 
a two-channel eikonal fit to all the high-energy pp and pp total cross section data above = 10 
GeV and to the elastic differential scattering cross section for pp at y/s = 31, 53 and 62 GeV, 
and for pp at y/s = 31, 53, 62, 546 and 1800 GeV. We use total cross section data sets compiled 
by the Particle Data Group [17] and differential cross section data sets from the references [18] 
(ISR), [19] (SppS) and [20] (Tevatron), see the review in Ref. [21]. The statistic and systematic 
errors of all scattering quantities have been added in quadrature. 

We introduce a parameter 7 to define the two diffractive eigenstates k of the eikonal model, 
such that their couplings to the Pomeron are 

(5 P , k {t) = (1± 7 )/3p(0, (29) 

where the proton wave function \p) = (|+) + |— ))/y / 2, see Ref. [14]. The CERN-ISR measure- 
ments of the excitations into particular channels (Nn, Nun, AK etc.) with M < 2.5 GeV [22] 
correspond to a cross section for low-mass diffraction of 7 

4°™ M ~ 2 mb at ^ = 31 GeV. (30) 

This value corresponds to 7 ~ 0.55, which we take in this analysis. 

We describe the pp and pp total and elastic cross section data in terms of Pomeron, and 
positive (j'2, 02) and negative (ou, p) signature Regge exchange. The positive signature secondary 
Reggeons, f'2 and 02, are taken to lie on an exchange-degenerate linear trajectory with intercept 
a+(0) and coupling 8 to the proton of parametric form 

(3 + (t) = (3 + (0) exp(r + t/2), (31) 

see (14). Similarly, the exchange-degenerate negative signature pair (uj,p) are described by 
parameters a_(0), /9_(0) and r_. The slopes, a' + , a'_, of the secondary Reggeon linear trajec- 
tories are fixed at 0.9 GeV -2 . On the other hand, the Pomeron coupling and trajectory are 
taken to have parametric forms given by (18) and (19), respectively. 

The values of the Regge parameters determined by the global fit are listed in Table 1. The 
fit has x 2 /DoF = 259/(325 — 11) = 0.83. The description of a to t and da/dt is displayed in Fig. 1, 
where the solid (dashed) curves correspond to pp (pp) scattering. The total cross section at 
LHC is predicted 9 to be a to t = 94.8 mb. 

7 Hcre, and in what follows, the value of ctsd accounts for dissociation of both colliding particles. 

8 If we assume the additive-quark-model relation, f3(f 2 ) — 3/3(a2), then the coupling (3(f 2 ) = \/9/10 (3 + . 

9 As expected, the inclusion of absorption reduces the value of the cross section predicted at the LHC energy. 
Indeed, if the high-mass diffractive absorptive effects were also included then the total cross section would be 
further reduced to about 90 mb [3, 23], see also [24, 4]. The present model also predicts an elastic pp differential 
cross section at the LHC energy with a diffractive dip at —t ~ 0.4GeV 2 . However, for this prediction to be 
reliable we must await the inclusion of the high-mass absorptive effects in the elastic analysis [23]. An earlier 
analysis [3], which included high-mass diffraction phenomenologically, showed that, at the LHC energy, the dip 
occurred at larger — t and that the elastic cross section was smooth in the region shown in Fig. 1. 
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Figure 1: Two-channel eikonal model description of total and differential cross section data. 
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/3p(0) 


2.26 
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0.02 


a,\ 


0.44 
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0.01 


a>2 


16.3 


± 


2.8 


a P (0) 


1.121 


± 


0.001 


OCp 


0.033 


± 


0.002 


P+(0) 


5.9 


± 


1.2 


r + 


0.5 


± 


1.2 


a+(Q) 


0.54 


± 


0.04 


/3_(0) 


2.4 


± 


0.9 


r_ 


3.1 


± 


6.0 


a_(0) 


0.57 


± 


0.09 



Table 1: The values of the Pomeron, and the positive and negative signature secondary Reggeon 
parameters obtained in the fit to the total and elastic differential cross section pp and pp data. 
GeV units are used; so, for example, the couplings /3(0) have units of GeV -1 . The couplings /3j(0) 
that we quote are actually Im^j/^O)), where the signature factor rji is given by (17), see the last 
paragraph of Section 2.1. The errors correspond to a 90% confidence level. 

4 High-mass diffraction: triple-Regge formalism 

The multichannel eikonal is unable to account for diffraction into high mass states. These 
processes, pp — ► pX with large Mx, are usually described in terms of a triple-Regge formalism, 
where 

M\ ~ (l-x L )s, (32) 

where Xl = 1 — £ is the momentum fraction of the ingoing proton carried by the outgoing 
proton. The approach is sketched in Fig. 2. Since we have large Mx we no longer have 




Figure 2: The triple-Regge description of high-mass diffractive dissociation, pp — > pX. In general, 
we may have vertices with three different Reggeons with coupling g^j. However, here it is sufficient 
to consider only contributions with i = i' . 
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b 2 _b± 
i 



Figure 3: A schematic diagram showing the notation of the impact parameters arising in the 
calculation of the screening corrections to the iij triple-Regge diagram. The conjugate momenta 

to 61, b 2 , b 3 are k t , q t , q t - If h = 0, then q t = q t . 



—t = q 2 . Rather, we must allow for non- vanishing t 

_ t = Qt , rn 2 p (l-x L y 

XL X L 

To the best of our knowledge, the screening corrections have not been explicitly included 
in the triple-Regge formalism. Therefore we present the formalism below, first using a single- 
channel eikonal and then generalising it to the two-channel case. Finally we treat the corrections 
to the 7T7rP diagram as a special case. 



t m i n . (33) 

XL 



4.1 Screening corrections in the triple-Regge formalism 

If we first neglect the screening correction, then the iij triple-Regge diagram of Fig. 2 gives 
the contribution 

Hup = (-srj ' (34) 

We use a simple exponential parametrisation of the triple-Regge vertices 

Quoit) = 9«M eMb'ujiq 2 + f~ k 2 t )), (35) 

where the momenta are defined in Fig. 3, and where q 2 = t m i n — q 2 / XL- 
Screening effects are best included by working in impact parameter space and using sup- 
pression factors of the form exp(— Q(b)). To determine the qt or t dependence we take the 
Fourier transforms with respect to the impact parameters specified in Fig. 3. We obtain 10 

M 2 da A f d 2 b 2 . , f d% «x . . f d 2 h 



dtdM 2 



A J ^ e ^ FM J ^ e ^F l{h ) J i^fo), (36) 
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Note that e ikt ' bl = 1 as k t = 0. 



where 



1 r I ' M 2 \~ a i l 

w-s^oi/^wb (38) 



2nPi(q t = 0) 

1 t f M 2 \~ a ^ 

■> 77 \ / d 2 k t l3Ak t ) — 

27T(3 j (k t = 0)J tHAtJ \s 

and where the ^-independent factors are collected in A 

/ s x 2a;(t min )-2 / M 2 \ aji0) ~ l 

a = mm^m (^) (-) . (39) 

These equations assume that we have a strong triple- Pomeron scenario, see (4). For the weak 
triple-Pomeron coupling ansatz, (2), we must include a factor q t in the integrand of the expres- 
sion (37) for Fj(6 2 ) when i,j = P, and also q t (= qt) in the analogous formula for Fj(6 3 ). 

To include the screening corrections, for a single-channel eikonal, we must include in the 
integrands on the right-hand side of (36) the factors 

exp(-Q(b 2 - 60/2) exp(-fi(6 3 - h)/2) = S(b 2 - h) S(b 3 - b ± ). (40) 

That is, we need to compute 



M 2 da 



dtdM 2 

where Id is given by 



nj 



= A/^F,(60IW| 2 , (41) 



7 d (&0 = j^e^F^Siih-h). (42) 



4.2 Generalisation to a two-channel eikonal 

Recall that in the two-channel eikonal, we take the Pomeron couplings to each diffract ive eigen 
component k to be 11 

p P , k (t) = (l± 7 )/3p(i). (43) 

On the other hand, for the secondary Reggeons, R, which are believed to dominantly couple to 
valence quarks, we take 

P R , k (t) = Pr®, (44) 

that is, the same /3r for each component. Thus for the PPP, RRP and ttttP couplings (with 
j = P) we make the replacement 

\h? - [(l + 7)|/ P | 2 + (l-7)|/ m | 2 ]/2 (45) 

11 Following [14], we assume the same structure, that is the same shape and size for each component, apart, 
of course, from the cross section. 
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in (41), whereas for the PPR and RRR (with j = R) we let 



\I p \ 2 + \I m \ 2 ]/2. 



(46) 



The subscripts p, m are used to denote the (1 ± 7) eigen components respectively. In addition 
for I p we replace S i: with % = P, in (42) by 



and for /„ 



S pP - - 

a _ 1 
iJmP — ~ 



;i+7)e 



-(l+ 7 ) 2 n/2 



+ (l- 7 )e 



-(l- 7 2 )n/2' 



(1 + 7 ) e -( 1 ^ 2 )^ 2 + (1 - 7 ) e -(i-7) 2 n/2 



(47) 



(48) 



These replacements of Si are for the PPP and PPR contributions. For the RRP, ttttP and 
contributions, with i = R or n, we have 



5* 



V - 2 

c _ 1 
3 m R — - 



( e -(l+ 7 ) 2 n/2 + e -(l-7 2 )^/2 



e -(l-7 2 )H/2 + e -(l-7) 2 H/2 



(49) 
(50) 



4.3 The screening of the ttttP contribution 

We treat the nnP contribution as a special case due to the presence of the n spin-flip amplitude. 
First, neglecting screening, we have 

MHo & N -t p.^H/MT^ 1 

dMAP = 4vr (T^ M)9mP \M*) [WJ ' ( ) 

where G 2 t1tN /Aix = 13.75 [25]. Since the pion is almost on-mass-shell, we expect g n7T p ~ |/?p(0), 
according to the additive quark model. In other words 

W//M0) ^ a(7rp)/a(pp). (52) 
The factor —t in (51) is given by 

. 2,2 2 , 2^- ~ X L . m S\ ~ x h) , , 

-t = Qt+ Q\\ = Qt+Qt + — » ( 53 ) 

where the qf = (eft ■ <?) 2 term is of proton spin-flip (sf) origin, and the remaining two terms 
are spin non-flip. Here, a is the Pauli matrix associated with the spin of the incoming proton. 
The non-flip interaction, is screened in the usual way, see (45), (49) and (50). However, the 
spin-flip amplitude is now a vector directed along q t . Thus the integrals in (37) and (42) should 
be written in vector form, with 

* (54) 



= 0) t-mV 
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and (/^ f (0)) 2 in (39) replaced by G 2 ^ n /Ati. This leads to 

F?(h) = b 2 rj(b 2 ) (55) 

and 

\lf\ 2 = \I*\ 2 + \Iy\ 2 (56) 

with x directed along q t . After the azimuthal integration, it means that the Bessel function 
Jo(qb), that arises in (37), should be replaced by J\(qb). 



5 Triple- Regge analysis of pp —> pX and pp —> pX data 

The main objective of this paper is to perform a triple- Regge analysis of the available d 2 a /dtd£ 
data for pp — > pX and pp — > pX (where £ = M|/s), allowing for screening (that is absorptive) 
effects. In this way we obtain a more reliable estimate of the bare triple-Pomeron coupling, 
gppp, which is a vital ingredient in the description of diffractive processes at high energy, see 
Ref. [12, 13, 26]. As can be seen from the triple- Regge formalism of the previous section, it is 
necessary to work in impact parameter space in order to include screening corrections. However, 
it is computationally time-consuming to evaluate the required Fourier transforms inside a \ 2 
fit. To facilitate the fit to the data, we therefore first prepare grids of the integral in (41) for 
a range of fixed values of the slopes b\^ of the triple-Regge vertices for each q t and M 2 data 
point. A MINUIT fit to the data is then performed using Chebyshev polynomial interpolation 
in b' Uj . 

In fact, we performed three triple-Regge fits to the d 2 a/dtd£, data with different assumptions 
for the form of the triple-Pomeron vertex. The first two assume either a strong or weak coupling 
triple-Pomeron vertex, that is they correspond to 

Sfp(0) and gY P (0)q 2 exp(-b'q 2 ) (57) 

respectively, see (4) and (2); recall that we use GeV units. The third fit considers a combination 
of the two couplings with both <?fp(0) and gfp(0) as free parameters. 

We fit to CERN-ISR 12 [27], FNAL fixed-target [28] and Tevatron [29] data for pp -> pX and 
pp — > pX. The differential distributions for the FNAL fixed target and Tevatron experiments 
can be found in Ref. [9] . The relative normalisations of the data sets are not very well known. 
In Ref. [9] it was claimed that the normalisation uncertainties of the FNAL data are about 
10%. We fix the normalisation of the ISR data and introduce normalisation factors for the 
FNAL data with a 10% error. However, we find that the Tevatron data prefer to be normalised 
up by the strong coupling fit, and down for the weak coupling fit, by factors of about 25%. We 
restrict this freedom and limit the normalisation factors to ±15% and ±10% for ^/s = 546 and 
1800 GeV data respectively, where the ± and — signs refer to the strong and weak coupling fits 

12 We choose a subset of the ISR data which is sufficient to fully describe their t and £ dependence. 
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respectively. The remaining normalisations coming from the fits are +10% for the yfs = 14, 20 
GeV data [28]. 

All three fits to the data prefer very small, or even negative, slopes, b' Uj , of the triple- Regge 
vertices. To avoid unphysical negative slopes, we impose the condition that b'^ > 0. In fact, the 
optimum x 2 for each fit has, within error bars, all b' U j = 0. Therefore, for the couplings quoted 
in Table 2, we have set all the slopes b' u ^ = 0. The only exceptions to this are the slopes of the 
PPP and PPR vertices for the weak coupling fit, which are found to be positive, as shown in 
Table 2. We fixed the nnP triple coupling at the additive quark model value, g nn p = |/3p(0). 
Actually the data prefer a smaller value, but with a large error. However the fixed value of g n7T p 
only enlarges x 2 by less than 0.5. The couplings of the Reggeons to the proton are taken from 
the "elastic" analysis of Section 3. The values of the remaining parameters, corresponding to 
the optimum triple- Regge fits, are given in Table 2. 





strong 


weak 


combination 


gip 
gfp 

UW 
Uppp 

gppR 

UW 
U PPR 

gpRp 

gRRR 


0.44 ± 0.05 

0.75 ± 0.10 

1.1 ± 0.3 
2.6 ± 1.0 


3.0 ± 1.2 
1.15 ± 0.3 
0.76 ± 0.15 
1.4 ± 1.7 
1.3 ± 0.5 
2.9 ± 1.4 


0.44 ±0.10 
0.2 ±0.5 

0.67 ± 0.16 

1.0 ± 0.5 
2.8 ± 1.5 


X 2 /DoF 


0.83 


1.40 


0.83 



Table 2: The values of the "bare" triple-Regge couplings guj(0) of (39), and slopes b' Uj of (37, 38), 
obtained in the three optimum fits to the cPa/dtdt; data. Recall that all the slopes b' Uj are set to 
zero, except for those of the PPP and PPR vertices in the weak coupling fit. 

We see from the \ 2 P er degree of freedom (DoF) that the data clearly prefer the strong, 
rather than the weak, triple-Pomeron coupling ansatz. This is also clear from the "weak±strong 
combination" fit, which is dominated by the strong coupling component, and is little changed 
from the pure strong coupling fit 13 . The preference for the strong triple-Pomeron coupling is 
also evident from Fig. 4. This compares the "strong" and "weak" descriptions of a sample of 
the d 2 a/dtdt; data that are fitted, including especially the high energy FNAL data at = 
546 and 1800 GeV. The FNAL data are most relevant since the triple-Pomeron contribution 
dominates for £ < 0.03. We see that the weak coupling description of the £ shape of these data 
is disfavoured. 

The t-dependence of the single diflractive cross section, d 2 a/dtd^, is shown, for example, in 
Fig. 5 for £ = 0.01, 0.1 at y/s = 1800 GeV, and in Fig. 6 for £ = 0.02, 0.06 at the much lower 

13 The normalisations of the data that are found in the "combination" fit are the same as those of the strong 
coupling fit. 
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Figure 4: The description of a sample of the d 2 a/dtd^ cross section data that are fitted using 
the strong (continuous curves) and weak (dashed curves) triple-Pomeron coupling ansatzes. (£ ~ 
M 2 1 s). Here, the curves corresponding to the (strong, weak) coupling fits of the FNAL data have 
been normalised (down, up) by 15% at y/s = 546 GeV and by 10% at y/s = 1800 GeV, to allow 
for the normalisations found for these data in the respective fits. 
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energy corresponding to s = 550 GeV 2 . We also show the individual triple- Regge contributions 
in these plots, together with the data available at these kinematic values. Note that, after 
accounting for the absorptive effects, the triple-Pomeron contribution does not vanish in the 
forward direction even in the weak coupling fit. Therefore to choose between the weak and 
strong coupling scenarios we have had to perform a full fit of the data. It is not enough just 
to study the low q t behaviour of the cross section. The first indication of the weak coupling 
characteristics is seen in Fig. 5 for £ ~ 0.01 and \t\ < 0.1 GeV 2 at the Tevatron energy. In this 
domain the data favour the strong coupling scenario. For larger values of £ the dip produced 
by the PPP term in the weak coupling case is filled in by the peak due to the irnP term. 
At the lower energy, corresponding to 550 GeV 2 , we see, from Fig. 6, that it is difficult to 
distinguish between the two scenarios. However, at the higher energy of the LHC the weak 
coupling scenario will reveal itself by a well pronounced dip at t ~ —0.02 GeV 2 after which 
d 2 a/dtd^ increases with \t\ up to \t\ ~ 0.15 GeV 2 , see Fig. 7. It is clear from Fig. 7 that 
single diffractive dissociation data obtained at the LHC for £ ~ 0.01 should be able to readily 
distinguish between the strong and weak triple-Pomeron coupling scenarios. 

In the remainder of this Section, we discuss only the favoured strong triple-Pomeron coupling 
fit, for which the triple-Regge couplings are the first set listed in Table 2. The values of these 
couplings are rather stable. They never go outside the quoted errors when we change the 
renormalisation parameters of the data, or allow g nwP to be a free parameter. The quality of 
the fit is excellent, with a minimum x 2 /DcF = 167/ (210 — 8) = 0.83. The 8 parameters are the 
4 couplings, gaj(0), and the 4 data renormalisations. The quality of the description can been 
seen from the plots of samples of the data that are fitted, Fig. 8 and Fig. 9, which also show 
the various triple-Regge contributions to the cross section. 

In comparison with the old triple-Regge analysis of Ref. [6], we now obtain a more than 
twice larger relative contribution of the PPR term. This is mainly due to the inclusion in 
our analysis of the higher energy Tevatron data. The inclusive cross sections measured in the 
interval of £ ~ 0.01 — 0.03 at the Tevatron energies turn out to be about twice smaller than 
that measured at the low CERN-ISR energies. 

The "strong coupling" parameters of Table 2 are much closer to the bare triple-Reggeon 
couplings than those coming from the old fits [5, 6, 7] which did not allow for the screening 
corrections. In particular the value 14 . 

93P = 9ppp ^ 0.2 (3 P (58) 

is consistent with the reasonable extrapolation of the perturbative BFKL Pomeron vertex to the 
low scale region [30] . However, these are still not the true bare vertices. For these sizeable values 
of guj the effect of "enhanced" screening, shown in Fig. 10, is not negligible. Moreover, it is not 
enough to just take one iteration, that is, to repeat the "elastic" analysis, but now including the 

1 Note that this result is in reasonable agreement with the parametrisation of the renormalised triple-Pomeron 
amplitude of Ref. [9] which leads to <?3p ~ 0.15 ftp 
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Figure 5: The four plots show the t-dependence of the d 2 a/dtd£ at £ = 0.01, 0.1 and ^/s = 1800 
GeV obtained in the strong and weak triple-Pomeron fits respectively. 
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Figure 6: The four plots show the t-dependence of the d 2 a/dtd^ at £ = 0.02, 0.06 and s = 
550 GeV 2 obtained in the strong and weak triple-Pomeron fits respectively, together with the data 
available at these kinematic values. 
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Figure 7: The continuous curves are the predictions for the t-dependence of the d 2 a/dtd^ at 
£ = 0.01, 0.1 and y/s = 14 TeV obtained in the strong triple-Pomeron fit. The disfavoured weak 
coupling predictions are shown by dashed curves. 



effects of the triple- Regge high-mass absorption, obtain new Regge-proton parameters and then 
repeat the triple- Regge analysis 15 . Rather, it is necessary to sum up the series of multi-Reggeon 
diagrams. We have considered such a model [23], that is an extension of an earlier model [3] 
but now including non-zero a' P as a free parameter, as well as more t-channel exchanges. The 
tuning of this model shows that it is possible to obtain a good description of the data provided 
that the triple-Pomeron coupling is a bit larger 

g 3P ~ 0.25 p P . (59) 



6 Inelastic J/ip photoproduction 

In Ref. [26] it was pointed out that the observation of the process 7P — * J/ip + Y at large values 
of My offers, in principle, an opportunity to determine the triple-Pomeron coupling where the 
screening corrections are smaller than in the pure hadronic reactions, see also [32]. Unfortu- 
nately, the M Y distribution has not been measured yet. However there exists a comparison 
of the HERA data for the "elastic" photoproduction process, ^p — > J/ip + p with the proton 

15 Note that an eikonal-type model for absorption will still violate s-channel unitarity [31]. At large b, on 
the periphery, where the opacity 0(6) is small, the effect of eikonal screening is not effective. For these (large 
I) partial waves the contribution coming from diffractive dissociation becomes larger than the total inelastic 
contribution. To satisfy unitarity more complicated (enhanced) multi-Pomcron diagrams must be included. 
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Figure 8: The description of the CERN-ISR pp — > pX cross section, d 2 a/dtd^, data [27] obtained 
in the strong triple-Pomeron coupling fit. (£ ~ M 2 /s). The individual triple-Regge contributions 
are also shown. 
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Figure 9: The description of the d 2 a/dtd^, measured in fixed-target and collider experiments at 
FNAL [28, 29, 9], obtained in the strong triple-Pomeron coupling fit. The individual triple-Regge 
contributions are also shown. 
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Figure 10: An example of an "enhanced" contribution to high-mass diffraction. 
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Figure 11: The process of proton dissociation in diffractive J ftp photoproduction, 7+p — > J/^ + Y, 
which is described by a diagram with a triple-Pomeron vertex in which the rescattering effects are 
small. The dotted line would mean the diagram became an enhanced diagram. This contribution is 
small. 



dissociation data. The ratio, at the photon-proton centre-of-mass energy W = 200 GeV and 
t = 0, is [33]- [37]. 

da(ip^ J/ip+pj/dt v ; 

where the "inelastic" cross section has been integrated over the mass region My < 30 GeV. 
Since only the Pomeron couples to charm quarks, the ratio r is described by the PPP and 
PPR contributions. The cross section for J/ip absorption is rather small. So we may neglect 
the screening factor, and obtain 



dM 2 f s \2«p- 2 



ir(3p(0) J M 2 \M 2 J 



9ppp 



M 2 



a P (0)-l 



+ 9PPR 



/3p(o) I so y 



(61) 



Note that the Pomeron trajectories ap(0) and ap(t) in (61), that is in the triple-Pomeron 
diagram in Fig. 11, are not the same. The lower Pomeron «p(0) in Fig. 11 is the usual 'soft' 
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Pomeron; whereas the upper ones, with ap(t), include DGLAP evolution from a low initial scale 
\i = hq up to a rather large scale \i ~ Mj^ at the J/ip production vertex. The summation of 
the double logarithms (a s ln(l/a;) ln(/i 2 //iQ)) n leads to a steeper x-dependence and hence to a 
larger effective intercept for the trajectory &p(t) of the upper 'hard' Pomeron. Thus to evaluate 
the ratio r of (61) we take dtp = 1.18, which corresponds to the W dependence observed in the 
HERA data [34]-[37]. Taking the other parameters from Table 1 and the strong coupling fit of 
Table 2, we obtain for the two terms in (61) 

r = rppp + rpp R = 0.12 + 0.06. (62) 

This is consistent with the HERA data, within the uncertainties. 

We also use the weak triple-Pomeron coupling fit of Table 2 to evaluate the ratio r at 
t = —0.2 GeV 2 . This is the lowest —t value for which r is measured. The experimental value 
is 16 [33] 

r = 0.4 ±0.1, for t = -0.2 GeV 2 . (63) 
When we evaluate the weak PPP coupling and PPR contributions to r, we find 

r = rppp + rpp R = 0.14 + 0.07, for t = -0.2 GeV 2 , (64) 

a value about a factor two smaller than that observed. So again the strong triple-Pomeron 
coupling is favoured over the weak. 

Originally, the strong triple-Pomeron coupling was supported by an analysis [26] of the 
HERA J/ip photoproduction data. The problem is that, at present, the My dependence of the 
cross section has not been measured; only the value integrated over an M y interval is given. 
The large value of the triple-Pomeron coupling that was extracted in [26] assumed that the 
whole cross section originates from Pomeron exchange, and neglected the secondary Reggeon 
contribution. Here we have included the secondary Reggeon contribution. However, if the M y 
distribution were to be observed, it would then be possible to separate these contributions 
using just the J/ip photoproduction data, and to give an independent check of the results of 
our triple- Regge analysis. Clearly a dedicated measurement of the process 7p — > J/ip + Y, 
including especially the M Y and t dependence, would be extremely informative. 



7 Conclusions 

We have performed a triple- Regge analysis of the available pp — > p + X and pp — > p + X data, 
in which we account for absorptive effects. Thus, the couplings of the triple-Reggeon vertices 
extracted in this analysis will be much closer to the couplings of the original bare vertices than 

16 The measured ratio at t = —0.2 GeV 2 is larger than that of (60) at t = 0, since the cross section of the 
elastic process, — ► J/ijj + p, has a larger t-slope than that for "fp — > J/ip + Y . 
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those obtained in the old analyses, where the screening corrections were implicitly included in 
the values of effective couplings. 

The triple-Pomeron vertex turns out to be rather large, g 3P ~ 0.2(3 P . This indicates that, 
in a more precise description of "soft" interactions, we must include more complicated 'en- 
hanced' diagrams and diagrams with the multi-Pomeron vertices (that is with a larger number 
of Pomerons in a vertex). 

Due to the inclusion of the high-energy Tevatron data in the analysis, we now obtain a twice 
larger relative contribution of the PPR term in comparison with the results of the old analyses. 

We considered two possible parametrisations of the small q t behaviour of the triple-Pomeron 
vertex: 

g 3P = constant and g 3P = fi^(0)g 2 exp(-6'g 2 ), (65) 

in GeV units. We found that the data prefer the first, the so-called strong Pomeron coupling 
scenario. In the second, the so-called weak coupling, case, we obtained a much poorer description 
of the high-energy Tevatron (CDF) data at small £ and low —t = 0.05 GeV 2 , where the triple- 
Pomeron contribution dominates. It will be important to measure the inclusive cross section 
for single diffraction, cPa/dtd!;, at the LHC to confirm this conclusion. This could be done 
when the forward detectors are operating at the LHC, even at moderate integrated luminosity 
[38]. Another possibility is to study in more detail J ftp diffractive production with dissociation 
of the target proton. That is, to measure the My and t dependence of the 7J9 — >■ J/ip + Y 
reaction. The HERA experiments have so far only presented the cross section integrated over 
My up to My = 30 GeV. However, already these data support the strong coupling solution. 

Finally, we note that our analysis of the data prefers zero slopes, corresponding to small 
size of the bare triple- Reggeon vertices. Together with the larger values of the PPP and PPR 
vertices, this gives hope that there is a smooth matching to the perturbative QCD treatment 
of the Pomeron. 
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